Huge metastability in high T c superconductors induced by parallel magnetic field 
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We present a study of the temperature-magnetic field phase diagram of homogeneous and inhomo- 
geneous superconductivity in the case of a quasi-two-dimensional superconductor with an extended 
saddle point in the energy dispersion under a parallel magnetic field. At low temperature, a huge 
metastability region appears, limited above by a steep superheating critical field (H. S h) and below 
by a strongly reentrant supercooling field (H sc ). We show that the Pauli limit (H p ) for the upper 
critical magnetic field is strongly enhanced due to the presence of the Van Hove singularity in the 
density of states. The formation of a non-uniform superconducting state is predicted to be very 
unlikely. 

PACS numbers: 74.60.Ec 
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The theoretical phase diagram of two-dimensional 
isotropic superconductors under parallel magnetic fields 
shows distinctive features such as a tricritical point where 
a high temperature second order phase transition gives 
way to a low temperature first order transition 0, Q and 
consequently to a region of metastability. Furthermore, 
in such systems, superconductivity may persist at mag- 
netic fields higher than the first order critical field if the 
formation of an inhomogeneous superconducting phase, 
i.e., the so-called Fulde-Ferrel (FF) phase, occurs. 

Superconducting copper-oxides have strong quasi-two- 
dimensional character, but the energy dispersion of the 
CuC>2 planes is clearly anisotropic and in particular, it 
contains extended saddle points as documented by an- 
gle resolved photoemission experiments. 5] Extended sad- 
dle points have also been found in the case of the two- 
dimensional Hubbard model by Assaad and Imada [|| 
using the Quantum Monte Carlo technique. Such sad- 
dle points lead to strong Van Hove singularities (VHS) 
in the density of states and recently, it has been sug- 
gested that such singularities could provide an expla- 
nation for the strange out-of-plane upper critical field of 
the High-T c superconductors. In this paper, we study 
the phase diagram of a two-dimensional superconductor 
with an extended saddle point in the energy dispersion 
under parallel magnetic field. 

In a strictly two-dimensional superconductor under an 
in-plane magnetic field, orbital frustration due to mag- 
netic field is not present and the upper critical field is 
determined by Zeeman pair breaking. Therefore, we con- 
sider the Hamiltonian 
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with = £fc — ah, = — /i and h = [ibH and where 
V, n, H, S and [Ib are respectively the attractive in- 
teraction constant, the chemical potential, the magnetic 
field, the area of the system and the Bohr magneton. 
The k sums in the interaction term follow the usual BCS 
restrictions. 

S-wave symmetry is assumed throughout the paper, 
for simplicity. The effects of d-wave symmetry in the 
case of a superconductor with constant density of states 
have been addressed by several authors. |9| It leads to 
visible modifications in the phase diagram for inhomoge- 
neous superconductivity with, in particular, the appear- 
ance of a low temperature kink in the phase boundary 
between the FF phase and the normal phase. However, 
the phase diagram for homogeneous superconductivity 
remains very similar to that of an s-wave superconduc- 
tor. 

Minimizing the free energy or, equivalently, fol- 

lowing a BCS mean field approach in the case of an 
homogeneous gap function, A = — VS^ 1 J2k( a - k -l ak ^) ' 
one obtains the two-dimensional gap equation 1 = 
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a £,k = V£k + ^ 2 an d where f(x) is the Fermi distri- 
bution function. Taking the limit A — *■ 0, one obtains 
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where Wrj is the usual frequency cutoff and t — kgT. 
From this equation, one extracts the temperature depen- 
dence of the critical field (H sc ) that induces the second 
order phase transition from the homogeneous supercon- 
ducting state to the normal state. Below a certain tem- 
perature (the tricritical point temperature), the phase 
transition becomes of first order and the field given by 
the previous equation becomes a supercooling field above 
which the normal state is a local minimum of the free en- 
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ergy. Note that this critical field has no dependence on 
the Fermi surface shape, only density of states depen- 
dence. 

In a two-dimensional system, a Van Hove singularity 
in the density of states results usually from the presence 
of a saddle point in the energy dispersion e(fe). In the 
case of a simple saddle point, e(q) ~ q 2 x — q 2 ., one has 
a logarithmic singularity in the density of states. In the 
case of an extended saddle point, e(q) ~ \q x \ n — \q y \ m , 
with higher powers and this leads to a power-law diver- 
gence in the density of states N(e) ~ (e — e v h)~ a with 

a = 1 — -. Besides two-dimensional systems, power- 

law divergences are present in one dimensional systems. 
In a one-dimensional system, a q 2 dispersion leads to a 
inverse square root divergence in the density of states. In 
the following, we will assume that the VHS is pinned at 
the Fermi level as observed in photoemission experiments 
in the copper-oxides Q and also indicated by theoretical 
studies of correlated two-dimensional models. [Tl| 

Assuming N(e) = N a \e — e v h\~ a , in the weak coupling 
limit, Eq. Elcan be rewritten as 
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with g = VN /2 and 
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The upper limit of the integral in Eq.[21was extended to 
infinity due to the fast decay of the integrand and there- 
fore, in the weak coupling limit, the supercooling field 
has no dependence on the frequency cutoff. The zero 
field critical temperature obtained from Eq.[3]is given by 
2^c0 ~ [g/( a ~ a 2 )] 1 /" which implies the well known en- 
hancement of critical temperature which motivated the 
Van Hove scenario of high T c cuprate superconductors. 
The zero temperature supercooling field (h sc ,o) is given 
by h sc , = (2 5 /a) 1 / Q .[l2 Therefore, h sc , ~ {l-a) l l a t c0 
and one has a much larger enhancement of t c o than that 
of h sc fi in the limit of very strong VHS, a — > 1. The 
reason behind the different enhancements is reflected in 
the fact that for fixed finite temperature and zero mag- 
netic field, the pairing susceptibility diverges in the limit 
a — *■ 1, but it docs not diverge for fixed finite magnetic 
field and zero temperature (note that tanh(x) is linear 
for small x). 

In Fig. ^ the temperature dependence of the reduced 
upper critical field (or supercooling field) obtained nu- 
merically from Eq. [2] is displayed for both a supercon- 
ductor with a VHS and a superconductor with constant 
density of states (BCS superconductor). The VHS ex- 
ponent a = 1/2 is obtained in the case of a quadratic 
one-dimensional energy dispersion or for example, in the 
case of an extended saddle point with quartic dispersion 
as found in the two-dimensional Hubbard model, jg] One 
observes in Fig^that the maximum second order critical 
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FIG. 1: The phase diagram of a paramagnetically limited two- 
dimensional superconductor with constant density of states 
(a — 0) and with a power-law divergence in the density of 
states (a — 1/2) pinned at the Fermi level. The reentrant 
behavior of the supercooling field H sc observed in BCS su- 
perconductors is strongly enhanced if a VHS is present in the 
DOS. 



field is not reached at zero temperature, but at an inter- 
mediate temperature. This reentrant behavior for the su- 
percooling field is known in BCS superconductors, Q and 
it has been recently observed in thin aluminum films. 
In the case of a VHS superconductor, this maximum is 
enhanced relatively to the zero temperature supercooling 
field and the reentrance becomes more pronounced. 

The strong reentrant behavior of h sc can be explained 
in the following way. At zero temperature, the Fermi 
surface splitting in the normal ground state due to the 
Zeeman coupling creates a no-pairing or blocking region 
around the zero field Fermi surface since, in order to con- 
tribute to the formation of a Cooper pair of momentum q, 
two states with opposite spins and of momenta k + q and 
— k must be either both empty or both occupied at zero 
temperature. In the case of homogeneous superconduc- 
tivity, Cooper pairs have zero total momentum, q = 0. 
At finite temperature, thermal excitations provide low 
energy pairing possibilities within this no-pairing region 
and therefore, the pairing susceptibility grows with tem- 
perature. In the presence of a VHS at the Fermi level, 
the number of these thermal excitations is much larger 
due to the strong splitting of the Fermi surface in the 
saddle point regions and therefore, this effect is more 
pronounced. 

Given a certain value of the superconducting gap, the 
free energy difference between the superconducting and 
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the normal state can be determined from 

F s (T,H)-F n {T 1 H)= f 
Jo 



(5) 
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The transition to the normal state with variation of field 
or temperature occurs when the zero gap local extreme 
of the free energy becomes the absolute minimum. If 
the finite gap local minimum of the superconducting 
phase converges to the normal state local extreme as the 
transition is approached, the transition is of second or- 
der, otherwise it is of first order and there is a region 
of metastability limited below by the supercooling field 
h sc and above by the superheating field h s h- The su- 
perheating field is the highest field for which there is 
still a finite gap solution of Eq. EI 13 It can be shown 
that below this field, there is a finite gap local mini- 
mum in the free energy. The first order critical field 
hi and the superheating field determined numerically 
from Eqs. and 03 are shown in Fig. ^ At zero tem- 
perature, the superheating field is just the zero temper- 
ature gap function, h s h,o — Ao = [/l^g] 1 /" ~ < c o, 
with A(a) — 27rcosec[(l — a)7r]P_ Q (0) where P u (x) is 
the Legendre function of the first kind. Consequently, 
h sc .o ~ (1 — a) 1 ^ Q /i s /i.o an d there is a huge metastability 
region. In Fig. ^ h is also apparent that the reduced 
temperature of the tricritical point is slightly larger in 
the case of the VHS superconductor. 

At zero temperature, the free energy difference is given 

by 
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The zero temperature critical field associated with the 
first order transition is usually denominated by Pauli 
limit or Chandrasekhar-Clogston limit [l4|. This field 
is easily extracted from the previous result, h p = {a(l — 

a)A(a)/4} 2^ A and therefore, this limit to supercon- 
ductivity is as strongly enhanced as the zero temper- 
ature superconducting gap due to the presence of the 
VHS. For a = 0, one recovers the well known result, 
h p = Ao/V^-Q In the limit a -> 1, h p -> A /2. Note 
that recent experiments indicate that the in-plane upper 
critical field in the copper-oxides exceeds considerably 
the BCS Pauli limit. pjj 

The previous results are independent of the shape of 
the Fermi surface. Now, we will consider finite q solutions 
of the gap equation, that is, we will search for a Fulde- 
Ferrel phase in the phase diagram. Following the BCS 
mean field approach in the case of an inhomogeneous gap 
function, A q = —VS~ l ^2 k (a-k-ql0>kt)i and taking the 
limit Aq — * 0, one obtains the following two-dimensional 
gap equation, 

V \p 1 ~ /(£p+q/2t) ~ /(£p-g/2j) 

s ^ 
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For a given temperature, the FF critical field is deter- 
mined by searching for the highest field solution of this 
equation for any value of q. 



It is well known that the Fulde-Ferrel state becomes en- 
hanced in the presence of Fermi surface nesting. [16] One 
might imagine that a VHS pinned at the Fermi level could 
provide a similar nesting effect in the Fulde-Ferrel state. 
However, one should be aware of a difference: while Fermi 
surface nesting does not enhance homogenous supercon- 
ductivity, in the case of VHS nesting it is considerably 
enhanced. It is therefore possible that the Fulde-Ferrel 
region of the phase diagram is narrower or that it even 
does not exist. We will show that at zero temperature 
and finite magnetic field, saddle points lead in fact to 
very poor Fermi surface nesting and indeed the Fulde- 
Ferrel region of the phase diagram is absent in the case 
of an extended saddle point. Furthermore, this phase re- 
mains absent even if one improves the nesting property 
by considering an isotropic Fermi surface with vanishing 
Fermi velocity. 

Let us consider the latter first. In order for an isotropic 
energy dispersion to have an VHS in the density of states, 
it must be of the form = — e v h — a-sign(g , )|g| where 
q = k — k v h- Again, we assume kp = k v h- The density 
of states for the above model is N(e) ~ a _1 ' b & _1 £^ _1 
and therefore the inhomogeneous gap equation can be 
rewritten as 
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with £± = [|£| 1/b sgn(£) ± l/2gcos6»] fc ± h where we have 
neglected a q 2 /k term since q <§; kp and where [x] b should 
be understood as sgn(x)|x| b . In order to keep the expres- 
sions simpler, we have considered a = 1. Note that in 
the one-dimensional case, the angle integration would be 
absent. 

In the weak coupling limit, the zero temperature criti- 
cal field, for q = q/(2h 1 / b ) < 1, can be obtained from the 
equation h c ,o(q) — [2gG(q)]& with 



JQ ~ J(l-qco h . 



2£ 



-i+i 



e)» * {£$ + qcos9) b + (£i - qcos9) b 

(9) ' 

Note that now the blocking region is no longer £ = h 
but £ = h{\ — <7Cos6>) h . G(x) has a maximum for 
x = 1 and therefore the maximum H c (q) is reached when 
q = q/{2h 1 / b ) = 1. This field is significantly enhanced by 
the VHS since it is proportional to V«, with a = 1 — i. 
However, there is only a weak enhancement in compar- 
ison with the q = supercooling field as observed in 
Fig. Q]in the case of a = 1/2. The FF phase bound- 
ary lies clearly below the first order critical field hi and 
therefore, this phase will not be observed. 

The same conclusion is reached if now one considers 
a energy dispersion with saddle points, = — e v h = 
Gdfel™ — l%| m ) where q = k— k vh , with the values of the 
momentum restricted to a small region around the saddle 
points k. vh and —'k vh by a cutoff. In this case, the zero 
gap pairing susceptibility (the integrand of Eq. [SJ has its 
maximum value for q = even for a simple quadratic 
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FIG. 2: The pairing susceptibility as a function of the renor- 
malized pair momentum q/(2h 1 ^ b ) for a BCS superconductor 
and a superconductor with an extended VHS at a small, but 
finite temperature, obtained numerically from Eq. |H| 



saddle point. In Fig. |3 plots of the zero gap pairing sus- 
ceptibility at zero temperature and for a fixed magnetic 
field are displayed for a BCS superconductor and for an 
extended saddle point with m = n = 4. The direction of 
momentum q for the latter is chosen to be along the x 
axis or along the diagonal of the Brillouin zone. For the 
saddle point, the maximum of the pairing susceptibility 
is for q = 0. In the case of a simple quadratic saddle 
point, the maximum remains at q — 0. If one "weakens" 
the saddle point by choosing exponents smaller than two, 
the maximum shifts to finite q, and rapidly becomes fixed 
at q = 1. One could therefore conclude that the Fulde- 
Ferrel state is absent in systems with energy dispersions 
containing saddle points close to the Fermi level. How- 



ever, one should be aware that the Fulde-Ferrel state is 
extremely sensitive to nesting properties of the Fermi sur- 
face and that we have only consider the saddle point con- 
tribution to the formation of this state and ignored the 
rest of the Fermi surface. In the case of homogeneous 
superconductivity, such procedure is justified since the 
saddle point contribution clearly dominates. In the case 
of Fulde-Ferrel superconductivity, if some portion of the 
Fermi surface is strongly nested, it may lead to a higher 
critical field than that of homogeneous superconductiv- 
ity. 

Concerning the experimental relevance of the results 
presented in this paper to the particular case of the 
cuprates, an obvious statement is that an huge metasta- 
bility region in the phase diagram is reflected by a low 
temperature hysteretic behavior which can be probed, for 
instance, by resistive critical field measurements |l7| or 
tunnelling measurements of density of states (as recently 
in thin Al films ^|). Unfortunately, in-plane critical 
fields of the high-T c superconductors are presently out- 
side the experimental magnetic field range. One can par- 
tially circumvent this difficulty by considering strongly 
overdoped or underdoped cuprates with lower critical 
fields. A study of hysteresis in these materials would 
provide a test of the validity of the extended Van Hove 
scenario for High T c cuprate superconductivity. 

In conclusion, we have shown that a two-dimensional 
superconductor with an extended saddle point in the en- 
ergy dispersion pinned at the Fermi level has at low tem- 
perature a large metastability region in the temperature- 
parallel magnetic field phase diagram and that the Pauli 
limit, Hp, for the upper critical magnetic field is strongly 
enhanced in this extended Van Hove scenario. Fulde- 
Ferrel superconductivity is absent from the phase dia- 
gram unless there are Fermi surface sections (away from 
the saddle point region) with very good nesting proper- 
ties. 
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